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[ We show that an exact non-perturbative quantization of continuum gauge the- 

2. ory is provided by the Faddeev-Popov formula in Landau gauge, 6{d ■ A) det[—d ■ 

^ ■ D{A)] exp[—SYM{A)], restricted to the region where the Faddeev-Popov operator is posi- 
Q-i' 

I I tive —d ■ D{A) > (Gribov region). Although there are Gribov copies inside this region, 

D . they have no influence on expectation-values. The starting point of the derivation is 



stochastic quantization which determines the Euclidean probability distribution P{A) by 
a method that is free of the Gribov critique. In the Landau-gauge limit the support of 
P{A) shrinks down to the Gribov region with Faddeev-Popov weight. The cut-off of the 
resulting functional integral on the boundary of the Gribov region does not change the 
form of the Dyson- Schwinger equations, because det[—d ■ D{A)] vanishes on the boundary, 
so there is no boundary contribution. However this cut-off does provide supplementary 
conditions that govern the choice of solution of the DS equations. In particular the "hori- 
zon condition" , though consistent with the perturbative renormalization group, puts QCD 
into a non-perturbative phase. The infrared asymptotic limit of the DS equations of 
QCD is obtained by neglecting the Yang-Mills action ^Sym- We sketch the extension to a 
BRST-invariant formulation. In the infrared asymptotic limit, the BRST-invariant action 
becomes BRST-exact, and deflnes a topological quantum field theory with an infinite mass 
gap. Confinement of quarks is discussed briefiy. 



1. Introduction 

Since the work of Gribov [1], a non-perturbative formulation of continuum gauge 
theory has appeared problematical due to the existence of Gribov copies. These are distinct 
but gauge-equivalent configurations A^"^^ = ^A^^'^ that both satisfy the gauge condition, 
d ■ A^^^ = d ■ A^'^^ = 0, where = g^^A^g + g~^d^g is a local gauge transformation. 
The difficulty arises when one wishes to quantize by gauge fixing namely by taking a 
single representative configuration on each gauge orbit. It has been proven that this 
cannot be done in a continuous way when space-time is compactified [2]. Geometrically 
this reflects the intricacy of gauge orbit space, the space of configurations A modulo local 
gauge transformations g. 

There is however an approach that by-passes the difficulties of Gribov copies by op- 
erating directly in A-space. This approach is stochastic quantization. For our purposes 
it is most conveniently expressed by the time-independent Fokker-Planck equation (given 
below) that determines the Euclidean probability distribution P{A). The geometric struc- 
ture of the equation assures that P{A) is correctly weighted. Although one cannot solve 
the Fokker-Planck equation exactly for finite values of the gauge parameter a, one can 
transform it into a system of Dyson- Schwinger (DS) equations for the correlation func- 
tions, that may be solved non-perturbatively, as has been done recently [3] . However these 
equations are more cumbersome than the DS equations in an action formalism. 

In sees. 2, 3, and 4, we find the exact solution of the time-independent Fokker-Planck 
equation in the Landau-gauge limit a ^ 0. The solution is remarkably simple. It is the 
familiar Faddeev-Popov weight, but restricted to the Gribov region fi,^ 

P{A) = N 6n{d ■ A) det[-d ■ D{A)] exp[-5yM(^)]. (1-1) 

The Gribov region is, by definition, the region in ^-space where A is transverse, and 
the Faddeev-Popov operator M{A) = —d • D{A) is positive, 

Q = {A: d-A = 0;-d- D{A) > }. (1.2) 

^ The Yang-Mills action is given by Sym{A) = (1/4) / d^x F^^ where F/^^ = 9^ A" - + 
gor'"'A''^Al, and the gauge-covariant derivative by [D^{A)u}]'' = d^Lu" + gof'^^Alio" . The 
Faddeev-Popov operator M{A) = —d ■ D{A) is symmetric when A is transverse, M{A) = 
-d ■ D(A) = -D(A) ■ d = Aft (A). Positivity of M{A) means all its non-trivial eigenvalues 
An (A) are positive. There is a trivial null eigenvalue with constant eigenvectors dfj,uj = 0, that are 
generators of global gauge transformations. In Appendix C we establish three simple properties 
of the Gribov region. 
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The first factor Sci{d • A) in (1.1) is the restriction of 6{d ■ A) to the region where M{A) 
is positive. Observables 0{A) are required to be gauge-invariant, 0{^A) = 0{A) and, 
by (1.1), expectation-values are calculated from 



where A^^ is the transverse part of A. Two comments are in order. 

(i) Gribov region Q vs fundamental modular region A. Formula (1.3) is paradoxical 
because the Gribov region fl is not free of Gribov copies [4] . The history of this formula 
is amusing. It was originally proposed by Gribov who conjectured in his seminal work [1] 
that there are no Gribov copies in fl. The same formula was also derived from stochastic 
quantization [5] by a method similar to the one presented in the present article (but using 
globally defined coordinates instead of coordinates defined only on a coordinate patch), 
and was interpreted to mean that the Gribov region fl is free of Gribov copies. However 
it was then proven [4], with details provided in [6], that there are Gribov copies inside fl. 
Moreover numerical studies [7], [8], [9], [10], and [11] revealed that in general there are 
many Gribov copies of a given configuration inside O. Consequently (1.3) was generally 
abandoned as an exact formula in favor of an integration over a region free of Gribov 
copies, known as the fundamental modular region A, 



The last formula is certainly correct and appears to contradict (1.3). It was subsequently 
argued nevertheless [12] that the functional integral (1.3) is in fact dominated by configura- 
tions on the common boundary of O and A. The derivation given in sees. 2, 3, and 4 shows 
that (1.3) is indeed correct. This is most fortunate because it is difficult to give an ex- 
plicit description of A. In Appendix A we examine concretely how the paradox is resolved. 
The lesson is that the normalized pobability distributions over A and O are equal in the 
sense that their moments of finite order n are equal. These are the correlation functions 
{A{xi)A{x2)---A{xn)) ■ This is possible in an infinite-dimensional space, where the proba- 
bility distribution may sit on a lower dimensional subspace such as a boundary or part of 
a boundary. This conclusion is consistent with numerical investigation of "Gribov noise" , 
namely the effect on measured quantities of taking different Gribov copies. Indeed for the 




N / dA'' 0{A'') ex-pi- Sym {A'')] det[-d ■ D{A'')], 



(1.3) 



{0{A))=N f DA'' det M{A'') 0{A'') exp[-SYM{A'')] 



(1.4) 
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gluon propagator in Landau gauge on reasonably large lattices, Gribov noise is quite small, 
of the same magnitude as the numerical accuracy [13], [14], [15]. The situation is quite 
different for a finite-dimensional integral, and the analogous problem for a finite lattice is 
also discussed in Appendix A. Formula (1.3) is also supported by a recent calculation in 
which the DS equation for the gluon propagator was derived from the time-independent 
Fokker-Planck equation at finite gauge parameter a. It was found to agree with the DS 
equation for the gluon propagator in Faddeev- Popov theory in the Landau gauge limit, 
a ^ 0, see particularly eqs. (9.4), (10.13), (10.14) and (10.17) of [3]. 

(ii) The form of the DS equations is unchanged by the cut-off on the boundary of Q. 
The DS equations are a set of equations for the correlation functions {A{xi)A{x2)---A{xn))- 
We shall derive them for the distribution (1.3) in sees. 5 and 6. They are compactly 
expressed as a single functional differential equation for the partition function or generating 
functional of correlation functions, 

Z{J) = N [ dA'' det[-d ■ D{A'')] exp[ -Sym{A'") + {J,A'') ]. (1.5) 

The functional DS equation for Z{J) follows from the identity, 

= nJjA'^ ^( det[-d ■ D{A'')] eM-SYM{A'') + iJ,A'')]), (1.6) 

which states that the integral of a derivative vanishes when there is no boundary con- 
tribution. There is in fact no boundary contribution, despite the cut-off on the bound- 
ary dfl, defined by the equation Ai(^*^) = 0, because the Faddeev-Popov determinent 
det[—d ■ D{A''^)] = Yin '^n(^^'^) vanishes on dCl. Thus the form of the DS equation is the 
same as if the integral were extended to infinity [16]. Again this is most fortunate because 
it means that implementing the restriction to the Gribov region causes no complication at 
all in the DS equations. 

Although the restriction to the interior of the Gribov horizon does not change the 
form of the DS equations, it docs provide supplementary conditions that govern the choice 
of solution. In fact the properties that result from the restriction to O, in particular 
the positivity of the weight P{A) and of the Faddeev-Popov operator M{A), dictate the 
natural choice of solution of the DS equation, that has been implemented previously, 
without necessarily invoking explicitly the cut-off at dfl, [17], [18], [19], [20], [21], [3], 
[22], [23], and reviewed in [24]. Another property is the horizon condition [25]. This 
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is an enhancement,^ compared to of the ghost propagator G{k) in the infrared, 

limk_,o[k'^G{k)]~^ = 0.^ In sec. 7 we show that the horizon condition is most conveniently- 
expressed as a formula for the ghost-propagator renormalization constant Z3. Although 
this formula flagrantly contradicts perturbation theory, it is nevertheless consistent with 
the perturbative renormalization group. The horizon condition puts QCD into a non- 
perturbative phase. 

In sec. 8 we deduce the asymptotic infrared limit of QCD by neglecting the terms in 
the DS equations that are subdominant in the infrared. It is found that the subdominant 
terms and only the subdominant terms come from the Yang-Mills action S'ym(^), so the 
infrared asymptotic limit of QCD is obtained by setting Sym{A) = 0. This is a continuum 
analog of the strong coupling limit of lattice gauge theory. The functional integral with 
exp[— S'ym(^)] replaced by 1 converges because it is cut off at the Gribov horizon. 

In Appendix B we outline the local BRST-invariant formulation of the present non- 
perturbative formulation. This assures that the Slavnov- Taylor identities hold at the non- 
perturbative level. In the infrared asymptotic limit, obtained by setting S'ym(^) = 0, 
the BRST-invariant action becomes BRST-exact, and defines a topological quantum field 
theory. As shown in sec. 9, this theory possess an infinite mass gap in the physical sector. 
In sec. 10 the extension to quarks is sketched out. 

The starting point of our derivation will be stochastic quantization of gauge fields. In 
the remainder of the Introduction we give a brief review of this subject so the reader may 
judge of the well-foundedness of this approach at the non-perturbative level. 

^ Entropy favors population near the boundary, in a configuration space with a high number 
N of dimensions, because of the volume element r^~^dr. The boundary dQ of the Gribov region 
n occurs where the lowest non-trivial eigenvalue of the Faddeev-Popov operator M{B) vanishes 
so, for typical configurations 5 on a large Euclidean volume V, M(B) has a very small eigenvalue. 
More precisely, compared to the Laplacian operator, M{B) has a high density per unit volume of 
eigenvalues p{X,B) at A = [25]. This enhances the ghost propagator G{x — y) = {Mxy{A)) in 
the infrared. 

^ The confinement criterion of Kugo and Ojima [26], [27], [28] yields the same condition in 
the Minkowskian theory. However for gauge-non-invariant quantities, the relation of the present 
approach, with a cut-off at the Euclidean Gribov horizon, to the Minkowskian theory remains to 
be clarified, perhaps along the lines of Appendix B. The relation of numerical gauge fixing by 
minimization in (Euclidean) lattice gauge theory to the Minkowskian theory is also not clear. 
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1.1. Review of stochastic quantization of gauge fields 

Historically, stochastic quantization originated [29] with the observation that the for- 
mal, unnormalizable Euclidean proabability distribution -Po(^) = ^exp[— S'ymI^)], with 
4- dimensional Euclidean Yang- Mills action 5'ym(^)) is the equilibrium distribution of the 
stochastic process defined by the equation, 

for the time-dependent probability distribution P{A,t). This equation is a continuum 
analog of the diffusion equation in the presence of the drift force Ki, 

dP d fdP 



(|^-K.P)=0. (1.8) 



dt dA^XdA 

that is known as the Fokker-Planck equation. If the drift force is conservative, Ki = 
— , then exp[— iS'ymI^)] is a time-independent solution. In Euclidean quantum field 
theory, t is an artificial 5th time that corresponds to the number of sweeps in a Monte- 
Carlo simulation, and that will be eliminated shortly. The same stochastic process may 
equivalently be represented by the Langevin equation 

where A^ = A'^{x,t) depends on the artificial 5th time. Here ry^ = r]'^{x,t) is Gaussian 
white noise defined by {v'^{x,t)) = and {7]l{x,t)T]'^{x,t)) = 25{x - y)Sf,J''^S{t - t')). If 
A'"exp[— jSymI^)] were a normalizable probability distribution — which it is not — every 
normalized solution to (1.7) would relax to it as equilibrium distribution. However the 
process defined by (1.7) or (1.9) does not provide a restoring force in gauge orbit direc- 
tions, so probability escapes to infinity along the gauge orbits, and as a result P{A,t) 
does not relax to a well-defined limiting distribution limt^co P{-^it) ^ A'"exp[— S'ym(^)] 
(although expectation-values of gauge-invariant observables formally do relax to an equi- 
librium value). 

A remedy is provided by the observation [30] that the Langevin equation may be 
modified by the addition of an infinitesimal gauge transformation, D'Jfv'^, 

dAl SS 



dt 



^ + + (1.10) 



where v'^ is at our disposal. This cannot alter the expectation-value of gauge-invariant 
quantities, for only a harmless infinitesimal gauge-transformation i^^gt,/i = -D^v has been 
introduced. In the language of the difi^usion equation, we may say that the additional drift 
force -K'gt,// is tangent to the gauge orbit. The modified Langevin equation is equivalent to 
the modified Fokker-Planck equation 

dt y ^^m«(x)Ua«(x) '''^^^^^ ^^^^^ 

We will choose v'^{x) to make D^v'^{x) globally restoring along gauge orbit direc- 
tions, so every normalized solution P{A, t) relaxes to a unique equilibrium distribution 
limt^^P{A,t) = P{A). 

Stochastic quantization in the time-dependent formulation has been developed by 
a number of authors who have expressed the solution as a functional integral [31], and 
demonstrated the renormalizability of this approach [32], [33]. A systematic development 
is presented in [34], [35], [36], [37], [38], [39], reviewed in [40], that includes the 4-and 
5-dimensional Dyson- Schwinger equation for the quantum effective action, an extension of 
the method to gravity, and gauge-invariant regularization by smoothing in the 5th time. 
Renormalizability has also been established by an elaboration of BRST techniques [41], [42]. 
Stochastic quantization may be and has been exactly simulated numerically including on 
rather large lattices, of volume (48)^, [43], [44], [45], [46], [47]. 

1.2. Time-independent stochastic quantization 

When the drift force is globally restoring, P{A) may be calculated directly without ref- 
erence to the artificial 5th time as the positive normalized solution of the time-independent 
Fokker-Planck equation 



HP 



f d^x— ^f--f^+K«p) =0 
J SA^^ix) V SA<^{x) ^ J 



(1.12) 

and Euclidean expectation values are calculated from (O) = J dA 0{A)P{A). We call H 
the "Fokker-Planck hamiltonian" . (It is not the quantum mechanical hamiltonian!). It has 
been proven directly [3] , without reference to the artificial time, that the expectation value 
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{0)y of a gauge-invariant observable A) = 0{A), is independent of v. Equation (1.12) 
determines a probability distribution P{A) directly in ^-space, that is correctly weighted 
at the non-perturbative level. The Gribov problem of globally correct gauge-fixing by 
identifying gauge orbits is by-passed. By contrast, in the Hamiltonian formulation of 
gauge theory, Gauss's law states that the wave functional ^{A) is gauge- invariant and is 
thus a functional defined on the space of gauge orbits [48] . 

To ensure that -K'gt,^ = D^v is globally restoring, we introduce a minimizing functional 
[49], [50], and [4], and choose -ftTgt,^ to be in the gauge-orbit direction of steepest descent. 
A convenient choice of minimizing functional"* is the Hilbert norm = J d'^x\A\'^. For 

an infinitesimal variation in the gauge-orbit direction = ^D^v^ we have 

5| I A| |2 = 2(A^, 5A^) = 2e(A^, D^v) = 2e(A^, d^v) = -2e{d^A^, v), (1.13) 

so steepest descent among gauge orbit directions of the minimizing functional is provided 
by V = a~^d • A with a > 0, and the time-independent Fokker-Planck equation is now 
specified to within a single gauge parameter, 

HP= I d^x—^ ( - — + KP) = 
J SA-(x)\ SA-(x) " J 



6A-{x) 



(1.14) 



(Symmetry and power-counting arguments also determine = a~^dxA'y^ = a~^d ■ A°'.) 

Having introduced the minimizing functional, we note that the Gribov region Q may 
be characterized as the set of relative minima^ with respect to local gauge transformations 
g{x) of the minimizing functional FA{g) = whereas the fundamental modular 



^ More generally, we may take for the minimizing function J d'^xA1^{x)af_ivA'^{x), where a^i, is a 
constant positive symmetric matrix. This defines a set of Lorentz-non-covariant but normahzable 
gauges that includes the Coulomb gauge as a limiting case [51]. To include different instanton 
sectors, one may choose as minimizing functional — y4„|p, where An is a fixed configuration of 
given instanton number. An alternative minimizing functional suitable for the Higgs phase was 
proposed in [42]. 

^ At any minimum, this functional is (a) stationary, and (b) the matrix of second derivatives is 
non-negative. These two conditions fix the properties that define the Gribov region: (a) transver- 
sality, d ■ A = 0, and (b) positivity of the Faddeev-Popov operator —D(A) ■ d. Property (a) 
follows from (1.13), which states that the first variation of the minimizing functional is 5||v4||^ = 
—2(u, d ■ A). Property (b) follows because the second variation is = —2{uj, d ■ D{A)ui). 
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region A may be characterized as the set of absolute minima. The set of absolute minima 
is free of Gribov copies, apart from the identification of gauge-equivalent points on the 
boundary dA, and may be identified with the gauge orbit space. In a lattice discretization 
the minimization problem is of spin-glass type, and one expects many nearly degenerate 
local minima on a typical gauge orbit, as is verified by numerical studies. Thus A is a 
proper subset of fi, A C fi, but A ^ O. 

1.3. Region of stable equilibrium of K^i 

The gauge transformation "force" Kgt is not conservative, and cannot be written, 
like the first term, as the gradient of some 4-dimensional gauge-fixing action, Kgt,ij. = 
a~^D^d ■ A'^{x) 7^ ^® cannot write the solution P{A) explicitly in general. 

However we shall solve (1.14) for P{A) exactly in the limit a ^ 0. In this limit P{A) gets 
concentrated in the region of stable equilibrium of the force i^gt,/^ = ci~^Dfj,d ■ A. 

Assertion: The region of stable equilibrium under the gauge transformation force 
Kgt,fj, = D^d ■ A is the Gribov region O. Proof: Transversality is a sufficient condition for 
equilibrium, because d ■ A = implies -K'gt,^ = 0. It is also necessary. For consider the 
fiow under this force, = D^d ■ A. We have d\\A\ \'^/dt = 2{A^, i^) = 2{A^, D^d ■ A) = 
2{An, d^d ■ A) = —2\\d • ^| P < 0, which is negative unless 9-^ = 0. We conclude that the 
region of equilibum under i^gt, which may be stable or unstable, is the set of transverse 
configurations. To find the region of stable equilibrium, observe that under this flow, we 
have §id-A = d-A = d-D{A) d-A. We linearize this equation to first order ind-A, which 
means taking d-D(A) d-D{A^') = -M(A^'), and we have f^d-A = -M(A^') d-A. Thus 
the equilibrium is stable when all eigenvalues of M{A^'^) are positive, and it is unstable 
otherwise. QED. 

2. A well-defined change of variable 

In order to solve the time-independent Fokker-Planck equation (1.14) in the limit 
a — > 0, we only need the solution for small a in a coordinate patch U in A-space that 
includes the Gribov region Q. In U, we make the change of variable A (B,g), defined 
by the gauge transformation, 

A^ = A^{B,g) = SB^ = g-^d^g + g-^B^g; with d-B = and M{B) > 0, (2.1) 
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where B e ^l. Local gauge transformations are parametrized by g{x) = ex-p[t°'9°'{x)] where, 
for each x, the 9°'{x) are coordinates for the SU(N) group. ^ The notation A = A{B,g) is 
understood to stand for A = A{B, 9), and we have B = A{B, 0). 

Gribov's critique of the Faddeev-Popov method is that this change of variable is not 
well-defined for all transverse B and g. We shall show however that it is well-defined in 
a coordinate patch U that includes f2. This is true, even though there are Gribov copies 
within f2, because the gauge orbits intersect f2 transversely. The coordinate patch U must 
be small enough in the ^-directions that the gauge transformations g{6) that relate these 
Gribov copies are not in U. 

To verify that the gauge-orbits intersect 0, transversely, it is sufficient to show that the 
change of variables (2.1) is invertible for infinitesimal angles 9"'{x) = e'^{x) for all B E Q. 
It follows that it is also invertible, and thus well-defined, on some finite cordinate patch U 
that includes O. 

To first order in e, the change of variable (2.1) is given by A^ = B^ + D^{B)e. 
The divergence of this equation reads d ■ A = d ■ D{B)e = —M{B)e, which shows that 
d ■ A depends linearly on e. Note that d ■ A is orthogonal to the trivial null space of 
M{B), consisting of constant functions, and we specify that e is also orthogonal to this 
null space. ^ Since S G O by assumption, M{B) is a strictly positive operator on the 
orthogonal space, and thus invertible, and we have e = -M-^{B) d ■ A. We solve for B 
in the form B^^ A^^ + D^^{B) M~'^{B) d ■ A. To zeroth order in e we have B = A = A^'^ , 
where A^^ = A^ — (9^)~^ d ■ A is the transverse part of A. This gives the inversion 
formulas B^ = Af, + D^{A^')M-^{A^') d ■ A and e ^ -M'^A^') d ■ A, vafid to first order 
in e or d ■ A. Thus for each A*"" G O, the change of variable (2.1) is invertible to first order 
in the small quantity d ■ A. QED 

Concerning the shape of the coordinate patch U, note that as the configuration B G 
Q approaches the boundary dQ of the Gribov region, the lowest non-trivial eigenvalue 
Xi{B) of the Faddeev-Popov operator M{B) approaches 0. Consequently the width in 
longitudinal or ^-directions of the coordinate patch U shrinks to zero as the boundary dO, 
is approached. We may picture W as as a very high-dimensional clam, shown in Fig. 1. 

^ Here and below we use the notation = f^A^^ and 5^ = f^B^j^. The are set of anti- 
hermitian traceless matrices that form the fundamental representation of the Lie algebra of SU(N), 
[t",t''] = f°'^H'', where the structure constants f"''"^ are completely anti-symmetric. 

^ The constant angles 8^9°' = parametrize global SU(N) transformations. These act within fl. 
However we may safely ignore them because they have finite volume that we normalize to unity. 
The spectrum of M(B) is discrete by quantization in a finite Euclidean volume. 
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3. Change of variable in Fokker-Planck equation 

To change variables in the Fokker-Planck equation, one takes over to functional vari- 
ables the standard formulas of differential geometry. The mechanics of the calculation 
are similar to the computation of the Coulomb hamiltonian by Christ and Lee [52], but 
there the change of variable was done globally whereas here it is done only in a coordinate 
patch. We freely go back and forth from continuum to discrete notation by the replace- 
ments A1^{x) ^ and (S^(a;), 0°'{x)) ^ u'^. In terms of A^, the Fokker-Planck equation 
reads, 

-HP ^ 4^5^' { 1^ - K^P] = 0, (3.1) 

dA' \dA^ / ' 

and expectation values are given by {F) = jY\.idA^ F{A) P{A). The coordinates A^ are 
Cartesian, but the coordinate transformation A = A{B, 6) = A{u) is non-linear, and the 
u — {B, 9) are curvilinear coordinates. In terms of these, the Fokker-Planck equation reads 



-HP ' ' 



0, (3.2) 



and expectation-values are given by {F) = /Ila'^'^"^/^ -^(^) -P(^)- The metric tensor 
is given by dA'dA' = du'^^^du^ = du'^Gc^pdul^ , with volume element VG = det 
The covariant and contravariant components of any Cartesian vector field Ki are given by 
Ka = l^i^i, and = ^K^. 

We now calculate these quantities explicitly in functional form. From A^^ = g~^Bg + 
g~^di^g, we obtain 

SA^ = g-' [sB^g + d^{6gg-') + [B, 6gg-'])g, (3.3) 

where 

io = dgg~' = ^ g-He^ = uj^dO^ = Tc^^d^'^ (3.4) 

is the Maurer-Cartan form. It satisfies dw = dgg~^ A dgg~^ = a; A a; or, in terms of 
components. 

We also have g~^t°'g = Raht^-, where the real orthogonal matrices Rah = Rab{0) = R^^ 

in the adjoint representation of the gauge group. From A^ = f^A^, and B^ = t"-B^, we 

obtain 

SA'', = Ral[K + Dt%^aSn], (3.6) 
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where is purely transverse, d^5B^ = 0, and D^f = D'^^{B) is the gauge-covariant 
derivative with the connection as argument. The last expression is the functional form 
of SA^ = ^^5u°^. It gives the functional operator that corresponds to and we have 
for the metric tensor, 

ds^ = J d^x 6A- 6A- = J d'x [5Bl + Dl^iu^Se-)] [5BI + Dl^iu^SO-)]. (3.7) 

To calculate y/G = det we start by writing the linear transformation (3.6) as the 
product of two transformations, SA'^^ = R~i^6C^, and 

sc; = 6B';, + D^;'iu^^50"). (3.8) 

The matrix Rat is orthogonal, so det i? = 1, and it is sufficient to calculate the deter minent 
of the linear transformation (3.8). We do this in two steps. We first transform from 5C^ 
to its transverse part (5(7)^"^'" = -P\^(5C)^, and its divergence SL"- = d^SC^, where = 
5xfj_ — d\{d'^)^^d^ is the projector onto transverse vector fields. This linear transformation 
is independent of the variables u = {B,9), so its determinent is a constant, and will be 
ignored. The linear transformation from SB and 69 to 6C^^ and 6L, is given by 

{5C)f^ = SBl + P^^Dl^^u^Jd-) 

where we have used the transversality of 5B^. This linear transformation is a triangular 
matrix, and its determinent is the product of the determinents of its diagonal submatrices. 
This gives 

VG = detI det[-df,D^{B)uj{e)] 

= det[-d^D^iB)] Detu{9) (3 ^0) 

= detM(S) JJdeta;(^(a;)), 

X 

which contains the Faddeev- Popov determinent det M(i?). It has been obtained by a 
purely local calculation at a fixed point A = ^B, without integrating globally over the 
gauge group. The volume element \/G is the product of det M(i?), that depends only 
on B, and the functional determinent Detcj(^) = *^®^'^(^(^))» ^^^^ depends only on 9. 
Here [det oj {9 {x))Yl^d9°'{x)] is the Haar measure of the SU(N) gauge group at x. It is 
common to write J Dg = J D9 Det(a;(^)). 
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We next find the inverse matrix |^ by solving for SB^ and SO'^. From (3.6) we obtain 

R,aSA- = [SB';, + D'l.^iuj'je-)]. (3.11) 
We take the divergence of this equation and use dfj^SB^ = to obtain 

d^{R,aSA';,)^d^D';,^{u;^je'^), (3.12) 

which gives the first inverse formula 

59- = J^[{d ■ D)-Yd^{RbaSA-^), (3.13) 

where J^{e) = {LJ-^)^{e). The Faddeev-Popov operator M{B) = -d ■ D{B) = -D{B) ■ d 
is symmetric and positive, so its inverse is well defined. To avoid a proliferation of indices, 
we write the last and similar equations in operator notation, 

69 = J{d-D)-^d-{R6A). (3.14) 

Inserting this into (3.11), we obtain the second inverse formula 

SBx = [5am - Dxid ■ D)-^d^]iRSA^). (3.15) 

One sees that 5Bx is transverse, d\5B\ = 0. The last two equations give the operators 
corresponding to the matrices From them we read off the continuum version of 

= M~ af^ namely, 

^=R{[5,x-d^{D.d)-'D,]^^-d,{D.d)-'j{e)^), (3.16) 



6A 



where R is the transpose of R. The {J-§Q)b = Jbi^)s§^ = i'^~^)b(^)w" angular 
momentum or Lie differential operators of the gauge group. They satisfy the Lie algebra 
commutation relations of the local gauge group 



s 



that follow from (3.5). 

We need the curvilinear components of the drift force = JCym,^^ + a~"^^^gt,/[i where 
-f^YM,^i(^) = ~ ^x^XniA) and -f^gt.^i = Df^d ■ A. We shall see that the one- form 

or covariant ^-component of Kym vanishes (because the action Sym{^B) = Sym{B) is 
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gauge invariant), while the tangent- vector or contravariant 5-component of Kg^ vanishes 
(because Kgt is tangent to the gauge orbit). Thus the Fokker-Planck equation (3.2) in 
curvihnear coordinates u= {B,9), reads HP = 0, where 

H = Hbb + Hbo + Hqg + Hqq, (3.18) 



(3.19) 



We use the continuum version of the formula KYM,iSA' = JB°' + K^^ JO'^ to 
obtain the one- form components of Kym- We have 

= j d'x D^FliB) [5Bl + D';^%aj^je-] (3.20) 
= j D^Fl^iB) 5B% 

by (3.6), where we have performed an integration by parts, and used {D^DxFxfj)"' = 
{l/2)gof"'^'^Fl^^F^^ — 0. Thus the one- form components of Kym are given by 

KYM,a = (4?l,a,4'ij = (I^A^^ (i^) , 0) • (3.21) 

We use the continuum version of Kg^^i-^ = K^fJ^g^ + K^^^^-^ to obtain the 
contravariant or tangent-vector components of -?^gt,/[i = -0^5 • A. We have 

^A^A = dxig-^Bxg + g'^dxg) = g'^ {dxidxgg'^) + [S, dxgg'^]) g 

= g-'Dx{B){dxgg-') g = g-^Dx{B){ujM'') 9, 

where we have used dxgg~^ = -^g~^dxO°' = Ua.dxO'^, and u is again the Maurer-Cartan 
form. In index and operator notation this reads 

^A^^ = Rah D\%uldxd'') ^ dxAx = R Dx{u;dx9), (3.23) 
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where D\ = D\{B). By the gauge transformation property of the gauge covariant deriva- 
tive D{A) = D{3B), this gives 

D^{A)d^Ay, = RD^{B) Dx{B){u;dxe), (3.24) 

By (3.16) we obtain 

/ ^lt,M^ = / MB) D^{B){u;dx9) 



(3.25) 

We perform an integration by parts and use Dn[5,j,i, — 0^(0 ■ d)~^Di,]-g^ = to obtain 

/ ^k.sJa = j pA(S)(a;a,^)]«[j(^)^]^. (3.26) 
Thus the tangent- vector components of Kg^ are given by 

^gt = = (0, J^,{e)[Dy^{B){ujdxet). (3.27) 



From (3.16) we obtain the Laplacian operator "^g^vG|^|^g^ in curvihnear 
coordinates, 



j d'x (^['^Am -Dx{d- D)-'d,] + ^J{e){d ■ D)-'d,) 
Putting all terms together, the explicit expressions for the terms in (3.19) are 



(3.28) 



X [5^. - d^{D ■ d)-^D,] —- - D,F,,{B) 



(3.29) 



X {d ■ D)-'[d, -d\D. d)-^D,] - D^F^^{B)] , 



(3.31) 



- - i?A[a'(^)9A^] )• 

a / 
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4. Solution in Landau-gauge limit 

We shall solve the Fokker-Planck equation HP = in the limit a ^ 0. In this limit 
the drift force in the gauge-orbit or ^-direction is dominant. This situation is reminiscent of 
the Born-Oppenheimer method in molecular physics. The 9 variables equilibrate rapidly, 
like the electron positions in a molecular wave function, and the dependence on the B 
variable is determined by an average over the 9 variable, like the nuclear variables. 

We expect that the solution gets concentrated close to 6' = 0. We rescale variable 
according to ^ = and find that Hbb is independent of a and unchanged, whereas 



and 



= ^ Det.(lv.e) / ^~(«^/^«) 



(4.2) 



are of leading order -^jj2 , while 



(4.3) 



is of leading order - . 

The Fokker-Planck hamiltonian has an expansion in a given hj H = oT^Hq + 
a-^/'^Hi + i?2 + 0(a^/2). We seek a solution of the form P = Pq + a^/^P^ + aP2 + 
which gives 

{a-^Ho + a-^/^Hi + H2 + ...) (Po + a'/^Pi + aP^ + ...) = 0. (4.4) 
To leading order we obtain 

-HoPo = j d^x-^[{d-D)-^ {-d^){D-d)-'^- D-c»e])Po = 0, (4.5) 
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or 



where D = D{B), M = M{B). The operator V = V{B) is defined hyV = M-^{-d'^)M-'^. 
It is symmetric and positive. 

The last equation is solved by a Gaussian in 0, 

Po(5,0) = QiB) N(detxy/^ exp[-(e,XG)/2] 

(4-7) 

= QiB) N{detxy/^ exp[-(^,X^)/(2a)], 

where {9,X6) = J d^x 6°'{x){X6)°-{x). Here X = X{B) is a symmetric operator to be 
determined, and N is fixed by 

j DO iV(detX)i/2 exp[-(^,X^)/(2a)] = 1. (4.8) 

The upper limit on the 9 integration actually finite, but this gives a correction of order 
exp(— 1/a) that we neglect. The solution (4.7) decreases rapidly as |^| increases away from 
0, as expected, with a Gaussian width \9\ ~ a^^"^. In the limit a — > 0, the support of the 
solution P{B,0) shrinks to ^ = 0, and is given by 

P{B,e)=5{e)Q{B). (4.9) 

We now check that (4.7) is actually the solution. Equation (4.6) yields two equations 
for X, 

(9, XV XQ) - (9, XM9) = 

(4.10) 

tT{VX - M) = 

that hold identically for all 9. The first equation yields 2XVX = XM + MX, or MY + 
YM = 2V for Y = X~^. Moreover when this equation is satisfied, the second equation is 
automatically satisfied. To solve for Y, we take matrix elements in the basis provided by 
the eigenfunctions of the Faddeev-Popov operator Mun — XnUn, and obtain 2{um, Vun) = 

{Xm + Xn){Um, YUn), OV 

{Um,X~^Un) = {Um,YUn) = 2{Xm + Xn)~^ {Um,VUn) 

(4.11) 

= 2 / dt {um,exp{-Mt) Vexp{-Mt) Un). 
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This gives 

poo 

X-^ =Y = 2 dt exp(-Mt) Fexp(-Mt) 

.00 (4.12) 

= 2 M"^ / dt exp(-Mt) {-d^)exp{-Mt) M'^ 
Jo 

and X = X{B) is indeed a positive operator, as is necessary for the normalizability of the 
Gaussian (4.7). 

The coefficient function Q{B) in (4.7) is left undetermined by the equation HqPq — 0. 
Since the leading term in the hamiltonian H — \Hq + ... leaves the solution indeterminate, 
we arc in the case of degenerate perturbation theory, and the lowest order solution is 
determined by a higher order perturbation. To obtain an equation for Q{B), we integrate 
the exact equation HP = over 0, 

J DQ Detio{a^/^Q) HP = 0, (4.13) 

where, we recall, H = Hbb + Hbo + Hqb + Hqq. This kills the Hq0 term that is of order 
^, for, by (4.3), it is the integral of an exact derivative, and thus vanishes identically, 
/ DQ Detcu(a^/20) HeeP = J DQ-^... = 0. For the same reason it kills the H^b term 
that is of order ^jy/j, / DQ DetCi;(a^/^0) HqbP = 0. It also kills the Hbo term that is of 
order because, by (4.1), the integral / DQ Beiuj{a}/'^Q) HseP is of the form 

/ D@ Detw(a^/20)J(a^/20)4Ti^ = - [ DQ Deta;(a^/20) F J(a^/^e)^l = 0, 
J oB J ()0 

(4.14) 

where the explicit form of F is not needed.^ The first equality holds by by the Lie 
group property that makes J(a^/^0)^ anti-hermitian with respect to Haar measure 
jDQ Deta;(ai/20). 

[It is easy to verify that the equation J DO Deta;(a-'^/^0) HboP = holds in the 
small-a limit. This is the same as the small angle approximation, and we have, to the order 



The fact that the integral on DQ surgically kills the Hqb and Hse terms is the pay-off for 
using the curvi-linear coordinates (B,9). In a previous calculation by the author [3], the time- 
independent Fokker-Planck equation was solved using Cartesian coordinates A^'^ and instead of 
{B,9). This gave an additional contribution, not surgically killed by the corresponding integration 
over DA^°, that was mistakenly neglected, and that was needed to cancel a spurious term, called 
K2, in the effective drift force. Fortunately K2 was neglected in [3], so what was thought to be 
an approximate formula there is in fact exact, and the calculation reported there is correct.] 
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required, g{e) = exp(^) = 1 + + (l/2)(t"^«)2. For the Maurer-Cartan form ^g'^ = 
f'u)^ we obtain, to the order required, uj^ = + (1/2)/«^^^t = 6""^ + «^/"7/307. The 
second term is an anti-symmetric matrix so for the Haar measure we get deta;(a^/^0) = 
1 + 0(a), and for the matrix jf , defined by J^u^^ = 6^, we get J^{a^/^Q) = + 



,1/2 



/a7/3eT + 0(a). This gives 



j DO Deta;(ai/2G)jf (a^/^Q) 



DO 



gal3 _|_ j^aifSQI \ 



) 



+ 0(a) 



The term in is an exact derivative because f°''^^ is anti-symmetric, and gives vanishing 
contribution. The leading term in F is of order so the remainder is of order a^/^ and 
vanishes in the small-a limit.] 

We conclude that in (4.13), the only surviving term is Hbb, given in (3.29). It is 
independent of a and 0, and (4.13) simplifies to 



Hbb Q = 0- 

From (3.29) we see that this equation is of the form 



(4.15) 



Q = 0. 



The left factor is orthogonal on v to longitudinal fields, so it may be written 



jtr 



...[S^..-d^{D-d)-'D,]P', 



5By 



+ 



Q = 0, 



where we have used the fact that functional differentiation with respect to a transverse 
field is ordinary functional differentiation with a transverse projector that comes from 



5Bl{x) 



ah 



(4.16) 



Thus the equation, Hbb Q{B) = 0, has the simple solution, 

Q{B)^Nexp[-SYM{B)]. 



(4.17) 



In continuum gauge theory, the Gribov region Q is convex, as shown in Appendix C, and 
therefore it is connected, so the normalization of the solution (4.17) is unique. We have 
obtained the solution in the coordinate patch U, in the limit a — > 0, 



P{B,9) = N6{d) e^-SYuiB)]. 



(4.18) 
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We express the solution P{B, 6) in terms of the original Cartesian coordinates A. 
The volume element is of course J dA. To first order in 9 we have A = B -\- D{B)9, and 
d- A = d- D(B)9, so 

6{9) = 6{d ■ A) det[-d ■ D{A)]. (4.19) 
Inside the coordinate patch U, the solution reads 

P(A) = N 6(8- A) det[-d ■ D(A)] exp[-S'YM(^)]. (4.20) 

Its support lies on 9 • ^ = 0, and it vanishes with det[—d ■ D{A)] on the boundary d^l of 
the Gribov region. We extend it to all of ^-space by stipulating that it vanishes outside U. 
For the diffusion equation with a drift force, the equilibrium distribution is unique [53] . 

5. Dyson-Schwinger equation for partition function 

To be of use, the non-perturbative Faddeev- Popov formula (1.3) must be supplemented 
with a prescription for how the functional integral, restricted to the Gribov region O, is to 
be evaluated non-perturbatively. An earlier approach [54] is to insert a ^-function 9{Xi{B)) 
that effects a cut-off at the Gribov horizon. The ^-function is given a suitable representation 
as an integral over auxiliary fields with a local effective action, and one integrates over all B 
without restriction and over the auxiliary fields. A far simpler approach [16] rests on the 
observation that the Gribov horizon dQ is a nodal surface of the integrand because the 
Faddeev-Popov determinent, det M{B) = ^n{B) vanishes with Ai(S), that is to say, 

on dVL. The DS equations, which are derived by a partial integration, do not pick up a 
boundary term, and would have the same form if the integral were extended to infinity. In 
this approach we never have to know where the Gribov horizon actually is. 

The partition function for the distribution (1.3) is given by 

Z{J) = N [ dB detM(S) exp[-SYM{B) + {J,B)], (5.1) 

where we have written B = A^^ , and (J, B) = J d^x J^{x)B^{x). Only the transverse part 
of J contributes, and we also take J to be identically transverse, J = J*"^. (The extension 
of the present non-perturbative approach with a cut-off at the Gribov horizon to an off- 
shell gauge condition with a local and BRST-invariant action is sketched in Appendix B.) 
The Faddeev-Popov determinent detM(S) vanishes on the boundary dil, so the identity 

^ = J^dB ^^i-^(^detM(S) exp[-SYM{B) + (J,S)]) (5.2) 
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holds, without any contribution from boundary terms even though the integral is cut-off at 
the Gribov horizon dQ. It is shown in Appendix C that the Gribov horizon surrounds the 
origin at a finite distance in all directions. 

To derive the functional DS equation for Z{J), we write det M[B) — exp[Tr In M(i?)], 
and define the total action 



S(S) = 5ym(S) -TrlnM(S), (5.3) 

so (5.2) reads 

^ = jjB (j^(a;)-||M)(^detM(S) exp[-5YM(S) + ( J, S)]) • (5.4) 

Although E(S) is not local in S, we shall derive the same DS equations as one gets from 
the usual local action of gluons and ghosts. We have 



= -[D^F^^{B)Y'^'^{x) - JgVM(^; B), (5.5) 



by (4.16), where "tr" means transverse part, [X^j\^^ = — 9^(9^) ^di,X^^ and the ghost 
current is given by 



J,^,^{x,B) _ - M [B) ) 



5Bl{x) V 5Bb^{x) 

_^^.^ S^OH-.gn)OA (M-^)--(.,.;i.)U.. (5.6) 



= -got"' j dS P;\{x-y) dx{M-'r{y,z;B) 



\z=y 



Here and below, derivatives act on the left argument of propagators. The identity (5.4) 
reads 

= [ dB U{x)+[D,F',^{B)f{x) + J^^Jx-,B)) 

^ ' (5.7) 

X (detM(S) exp[-5YM(5) + (J,S)]), 

and yields the functional DS equation for the partition function Z{J), 

-iDxFl[^)nx) Z{J) = [ J) + J» ] Z{J\ (5.8) 
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where DxF^^( jj) is a cubic polynomial in jj. The quantum ghost current in the presence 
of the source J is, by (5.6), 



(5.9) 



Here we have introduced the ghost propagator in presence of the source J, 

g^^ix, y;J) = { (M-i)-(a;, y; B) )j, (5.10) 
where {0)j denotes the mean value of 0{B) in the presence of the source J, 

{0)j ^ Z-\J) N [ dB detM{B)0{B) exp[-5YM(S) + ( J, 5)]. (5.11) 
Jo. 

To obtain a closed system of equations, we need a DS equation for the ghost propagator 
Q°'^{x,y; J). It contains a term proportional to X[^{B), so we must avoid integrating by 
parts on B or introducing ghost sources. (But see Appendix B.) Fortunately the functional 

DS equation for Q'^^{x, y. J) follows from the trivial identity / = M{B) M~^{B), that we 
average with P{B) exp[(J, B)], 



5{x - y)5'^' Z{J) = [ DB M'^^B) {M-')ll{B) P{B) exp[(J,S)] 

JO. 

= ^""{jj) j^DB {M-')%{B) P{B) exp[(J,i?)], 



(5.12) 



where M«'=(^) = - g^p^'^-A^d^,. Here P{B) = det M(S) exp[-5YM(S)] is the 

probability distribution, although the form of the DS equation for the ghost propgator is 
independent of P{B). This gives the DS equation for the ghost propagator 

M«^(^) [g''\x,y;J) Z{J) ]^5{x-y)5'^' Z{J). (5.13) 

Equations (5.8) and (5.13) and formula (5.9) provide a complete system of functional DS 
equations for the partition function Z{J) and the ghost propagator g^^{x,y; J). 
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6. Functional DS equation for gluon and ghost propagators 

We change variable from Z[J) — expW{J) to the "free energy" W{J). For the ghost 
propagator we obtain 

'SW S 



M 



ac 



We again change variables by Legendre transformation from the free energy W{J) to the 
quantum effective action 

r(Sei) = J,Sei,, - W{J), (6.2) 

where the new variable B^^ ^{x) is defined by 

It is identically transverse, -Bci,^ = B^l ^, and takes values in Q because Bc\{J) = {B)j is 
an average with a positive probability, N detM(S) exp(S, J), over the convex region Q. 
Inversion of B^i = -Bci(J) to obtain J = J{Bc\) is possible because the gluon propagator 
in the presence of the source J, 

V,y{j) ^ ( (S. - {B^)j) {By - {By)j) )j = = , (6.4) 

is a positive matrix. The gluon propagator is expressed in terms of the Legendre- 
transformed variables B and r{B) by 

= («-^) 

Here and below, we write B instead of B^. The gluon propagator and its inverse are 
identically transverse, d\D\^{x^y^B) = 0. 

Under the Legendre transformation, derivatives transform according to 

as one sees from (6.4). In terms of the Legendre transformed variables, the DS equation 
(6.1) for the ghost propagator reads 

5 



5{x-y)5-' = M'^'^[B + V—) G''\x,y;B) 
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M-%B)g^\x,y;B)-gor'''' / dz V^^^x, z; B)j^d^g-\x,y; B), 

(6.7) 



where Q^^{x,y; B) = Q^^{x,y; J{B)) is the ghost propagator expressed in terms of the 
source B. Finally, instead of g''^{x,y;B), we take as new unknown variable the inverse 
ghost propagator T^f^{x,y; B) defined by 

rgh,.,(5) = 6;-i(5) ^ J dyT'$,{x,y;B)g''%y,z;B) = 5{x-z)5--. (6.8) 



We substitute 



^ -Q.y{B) = -gUB)^^^^Qvy{B). (6.9) 



into the previous DS equation, and multiply on the right by the matrix T^\^^y^ to obtain 
the functional DS equation for the inverse ghost propagator 

r^^(a;,y;5) =M-\B)5{x - y) 

, f . f 6Ti'{u,y;B) (6-10) 

+ gor"' J dzdu V^lix, z; B)d^g'^{x, u; B)^^^^, 

that is represented diagrammatically in Fig. 2. Here — fB^{z) — complete ghost- 

ghost-gluon vertex in the presence of the source. 

We make the same changes of variable in the functional DS equation (5.8) for Z{J) — 
exp[VF(J)]. We evaluate^ 



5_ 

= exp[W^( J)] [dx (dxBl - d^Bl + gor'''{B\Bl + V\l{x, x, B)) 



^^^^"'^(^j)^^^ exp[M^(J)] 



r (6.11) 

+ gor'^ [bI + —^g,r'^^{BiBl + Vliix, X, B)^ 



= exp[W^(J)] ^xF^,{B) + Jj„.g,,^(a;; B)\ , 

where B^ = ^^j"^^ and V^i,{x^ y;B) = ^j"!^) • The quantum gluon current in the presence 
of the source B is defined by 



X / dydzdw V'/^{x, y; B)V'll{x, z; B)Vf^{x, w; B) Vi^f^iy, z, w; B 



2 fabc fcde 
tr 



(6.12) 



^ In this section we write VJi''(A) = dij,5°''^+gof°'^'^A''^ for the gauge-covariant derivative instead 
of D^°{A) to avoid confusion with the gluon propagator T). 
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Here 

rg^(y, z, w; B) = (6.13) 

is the complete triple-gluon vertex in the presence of the source B. This gives the functional 
DS equation for V{B), 

= - '^xFx,{B){x) - jJu.gh,M(^; B) - Jl.,,,,{x; B), (6.14) 

where, by (5.9), 

^qu.gh,p(^) B) = i^qu.gh,/i(^) J{B)) 

u f . (6-15) 

= -gor'" J d% P;\{x - y) dxG'^^iy, z; S)U=, 

is the quantum ghost current in the presence of the source B. 

A more explicit form of this equation, is obtained by differentiating with respect to 
B^{u), which yields a functional DS equation for the inverse gluon propagator, 

SbZTbUu) < - + - ^«°^""^-)"*(^ - ") (6.16) 

+ (ghost loop) + (1 gluon loop) + (tadpole) + (2 gluon loops) 

where 

(ghost loop) = -gor''[ J dydz d^g"'ix,y; B) g'^%x,z;B) (6.17) 

(1 gluon loop) = j dydz (gof^'' (Sx^S,,^ + S^xS^^ - '^SxJnfi) 

+ ^xT^U^.y^B) Vli{x,z;B))'\%%{y,z,u;B), 

(tadpole) = glr'^'f^'' ( {5x^.6^. + 5^x5^^. - 25a.5«m) ' K^i^^ B) (6.19) 
where superscript "tr" means projection onto transverse parts on (x, k) and {u,t). The 
complete ghost-ghost-gluon vertex in the presence of the source B, — sB^lu) — ' reappears 
in (6.17), and the complete triple-gluon vertex r^^^(j/, z, u; B) in the presence of the source 
B is defined in (6.13). We do not write out explicitly the two-loop term, but all terms are 
expressed graphically in Fig. 3. 

The pair of equations (6.10) and (6.16) are a complete system of functional DS 
equations for the quantum effective action r(S), and for the inverse ghost propagator 
T'^{x,y\B). These functional equations are converted to equations for the coefficient 
functions by differentiating an arbitrary number of times with respect to B, and then 
setting 5 = 0. 



(6.18) 
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7. Horizon condition and renormalization 



Solutions are subject to the supplementary conditions that both the gluon and ghost 
inverse propagators j^a(^-)^gi (-^■) cind T'^^{xjy; B) be positive matrices. Another supple- 
mentary condition results from the fact, discussed in Appendix A, that in a space of high- 
dimension, entropy favors a high concentration of population very near the boundary dfl of 
the bounded region Q. The boundary occurs where the lowest non-trivial eigenvalue of the 
Faddeev-Popov operator M{B) vanishes. Thus, for typical configurations S, the positive 
operator M{B) has a very small eigenvalue and, in fact, it has a high density of eigenvalues 
p(A, -B) at A = 0, per unit Euclidean volume V, as compared to the Laplacian operator [25]. 
This makes the ghost propagator, G{x — y)d°'^ = ((M~^)"^(x, y; -B)), long range, so in mo- 
mentum space it is enhanced at p = compared to limp^Q[p'^G{p)]~^ = 0, [1], [25]. 
This property will provide a non-perturbative formula for the ghost-propagator renormal- 
ization constant Z3 that moreover is consistent with the perturbative renormalization- 
group. 

The gluon and ghost propagators, with source B = 0, are given in momentum space 

by 

D^^{x) = (27r)-4 J d^k b^u{k) exp(z/e • x) 
G{x) = (27r)-^ J G{p) exp(zp • x), 
and the ghost-gluon vertex by 



(7.1) 



/ T,ix-y,y-z) = ^^^j-^\s=o 



(7.2) 



The DS equation for the ghost propagator G(p), obtained from (6.10) by setting S = 0, 
reads 

G-\p) =p^- Ngo ip^ (27r)-^ j d^k b^,{k) G{p - k) V,{p - k,p). (7.3) 

All quantities are unrenormalized, and we have used J°^^'^ycde _ jsj-^ae SU(N). 

Factorization of the external ghost momentum is a well-known special property of the 
Landau gauge that makes it less divergent than other gauges. To make it explicit, we 
note that the ghost-ghost-gluon vertex {p, q) is a function of two linearly independent 
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4- vectors. It is also transverse, {p — q)ijTij,{p, q) = 0, because the transversality condition 
is imposed on-shell, so it may be written 

f 9) = -igo P';M p. V{p\ q^), (7.4) 

where k = q — p. The scalar vertex function, V{p^ .k"^ ,q^) is symmetric V {p"^ , k'^ , q'^) = 
k'^,p^) in consequence of the symmetry Q'^'^{x, z\ B) = Q'^°'{z, x; B). The DS equation 
for the ghost propagator reads, 

G-\p) =p^- Ngl p^p, (27r)-4 j d^k D^,{k) G{p - k) V{{p - k)\ k\ p"), (7.5) 

where the factorization of the two external ghost momenta p^ and Pi, is now explicit. 

This equation is divergent and must be renormalized. In perturbative renormalization 
theory, quantities renormalize according to 

D^^^Z^Dr,^,- G = Z^Gr- V^Z^^Vr; go = Z^{ZsZl/y^gR, (7.6) 

and in Landau gauge the additional special property 

Zi = l; V = Vr; go = {ZzZ]!\^gR (7.7) 

holds. In terms of renormalized quantities, the DS equation for the ghost propagator reads, 

Gr\p) = p'Za - Ngl p^p, (271)-^ J d^k DR,^,{k) Gr{p - k) Vr{{p - kf, k\ p"). 

(7.8) 

To avoid infrared difficulties, the ghost propagator is usually renormalized at some 
finite renormalization mass jj.. However the horizon condition, limp2^ob^^(^')]~^ — ^i 
allows us to renormalize at p = 0. It tells us that in the DS equation (7.8), the first 
term, p^Z^-, must be cancelled by the term of order p^ in the second term. This gives a 
renormalization condition at p = 0, in the form of an equation for Z3, 

Z3 = Ngl (27r)-^ / d'^k p^,p.DR,^,Ak) GR{k) VR{k^ 0), (7.9) 

^|fc|<A 

where A is an ultraviolet cut-off. We have set p = in the integrand, and the integral is 
independent of the direction p. This statement of the horizon condition shows that it is 
flagrantly non-perturbative because, in perturbation theory, the left hand side is of order 1, 
but the right hand side is of leading order gR. 
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The last equation gives the renormahzation-group flow, 

= Ng\ {p^p./p") (27r)-^ ^^(A) Gr{K) Fi^(A^ A^0) 

= Ngl{A7:)-\Z/2) A^ ^«(A) Gr{K) Vr{A\K\Q). 

As a check, we note that if we take the tree values Dr[A) = Gr[A) = 1/A^, and 
V{p'^, /c^, g^) = 1, we obtain 

^^-i^^)-'m)N9'o + 0{g',). (7.11) 

The term of order is scheme-independent, and agrees with the standard one- loop expres- 
sion in Landau gauge. Thus the horizon condition provides a normalization condition for 
the ghost propagator at p = that is in flagrant disagreement with perturbation theory, 
but nevertheless satisfies the perturbative renormalization-group flow equation. 

We substitute (7.9) into the DS equation (7.8) for the ghost propagator, and obtain 

Gr\p) = Ngl p^p, (27r)-4 j d^k DR,^,Ak) 

X \GR{k) VR{k^e,0)-GR{p-k) VR{{p-kf,k^p% 

This is a finite, renormalized DS equation for the ghost propagator. It is invariant under the 
renormalization group in the sense that it is form-invariant under the transformation (7.6) 
and (7.7) of perturbative renormalization theory in Landau gauge. This equation, from 
which the tree term k"^ has been eliminated by the horizon condition, gives the ghost prop- 
agator an infrared anomalous dimension qq, so it behaves likes G{k) ~ {/j,'^)"''^ / {k'^Y'^"''^ 
in the infrared. This puts QCD into a non-perturbative phase. 



8. Exact infrared asymptotic limit of QCD 

Recent solutions of the truncated coupled DS equations for the gluon and ghost prop- 
agators yield ghost propagators that are enhanced in the infrared, and gluon propagators 
that are infrared suppressed [17], [18], [19], [20], [21], [22], [23], and [3]. Typical values for 
the infrared asymptotic form of the gluon and ghost propagators [20] and [21] are, 

G^ik) = //2'^-/(/c2)l+«« ^ (/x2)0-595/(;^2)1.595^ ^^'^^ 



27 



ac = (93 — -\/T201)/98 ^ 0.595, an = —'^clgi where ao and ao are the infrared critical 
exponents of the ghost and gluon. The gluon propagator D{k) is so strongly suppressed 
at A; = that it vanishes -D(O) = 0. With D{x — y) = {A{x)A{y))^ this corresponds to 
suppression of the low-frequency modes of A{x) in the functional integral. The actual 
values of the infrared critical exponents do not depend too strongly on the truncation 
scheme [20]. The salient infrared features are easily understood. The cut-off of the func- 
tional integral at the Gribov horizon is implemented in the DS equations by the horizon 
condition. It states that the ghost propagator G{k) is enhanced in the infrared or, equiva- 
lently, that the infrared critical exponent of the ghost positive, ac > 0. The DS equations 
yield ao = —2ao-, so enhancement of the ghost causes suppression of the gluon in the 
infrared. This is the expression in the DS equations of the proximity of the Gribov horizon 
in infrared directions. 

The results of calculation with the DS equations are in at least qualitative agreement 
with numerical evaluations of gluon and ghost propagators [55], [56], [13], [57], [58], [59], 
which, on sufficiently large lattices, yield a gluon propagator D{k) that turns over and 
decreases as k decreases [60], [61], [62], [63], [64], [65], [66], [67] (possibly extrapolating 
to i^(0) = at infinite lattice volume), with a turn-over point kmax that scales like a 
physical mass [68] . The only explanation for this counter-intuitive turn-over is the strong 
suppression of infrared components by the proximity of the Gribov horizon in infrared 
directions. The agreement of DS and numerical calculations gives us confidence that we 
have a reliable picture of the gluon and ghost propagators including, in particular, in the 
infrared region. 

One may use the above expressions for the asymptotic propagators to estimate the 
convergence and magnitude of the various terms on the right hand side of the DS equations, 
simply by counting powers of momentum. The dominant terms in the infrared region are 
the ones that contain the most ghost propators G{k) in the loop integrals. The infrared 
limit of the truncated DS equations are found to have the following remarkable properties: 

(i) The infrared limit of the DS equations decouples from the degrees of freedom as- 
sociated with finite momentum and is free of ultraviolet divergences. Technically, what 
is found is that when the external momenta kf, are small compared to Aqcd, then the 
internal loop momenta ki scale like the /ce, and the contribution when the ki are large 
compared to ke may be neglected. As a result, when the ke are small, one may replace the 
propagators and vertices in internal loops by their infrared asymptotic forms G^{k) and 
D^{k) etc. The loop contributions that are dominant in the infrared are the ones that 
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contain the most ghost propagators G{k).^^ The asymptotic infrared hmit of the DS equa- 



tions is highly convergent in the ultraviolet because G^^{k) is strongly suppressed there. 
In fact the asymptotic gluon equation, given below, is finite without renormalization, and 
the asymptotic ghost equation is finite with the renormalization (7.12). We conclude that 
the DS equations possess an infrared asymptotic limit that is well-defined, and decoupled 
from propagators and vertices at finite momentum. 

(ii) The terms that are dominant in the infrared limit come from the action 
— TrlnM(S), whereas the subdominant terms come from Yang-Mills action Sym{B). It 
is instructive to classify terms that are dominant or subdominant on the right hand side 
of the DS equations according as they originate with the action, —Tr In M(i?), or with 
the Yang-Mills action, Sym{B). Because the ghost propagator is enhanced in the infrared 
while the gluon propoagator is suppressed, one finds that all subdominant terms and only 
the subdominant terms disappear if one sets Sym{B) = in the derivation of the DS 
equations given in sees. 5 and 6. 

Because the solutions of the truncated DS equations are consistent with numerical 
evaluations of the gluon and ghost propagators, the effects of truncation should not be too 
drastic. We therefore expect that properties (i) and (ii) of the truncated DS equations hold 
also for the solutions of the exact, untruncated, DS equations, that is, that there exists an 
exact infrared asymptotic limit of the DS equations that is obtained by setting S'ym(-B) = 0. 
This implies that the cut-off at the Gribov horizon suffices to make the functional integral 
over B converge, even though exp[— 5'ym(-B)] is replaced by 1. 

We now write in functional form the exact infrared asymptotic DS equations (with- 
out truncation!) that are obtained by setting Sym{B) = 0. We designate the generat- 
ing functionals where the coefficient functions are given their asymptotic forms by T{B) 

For the ghost-propagator equation (7.3) or (7.12), both terms on the right-hand side are dom- 
inant, and both originate from the action — TrlnM(fi). The gluon-propagator equation (6.16), 
with source B — 0, reads 



The tree term, of order k^, is subdominanant in the infrared compared to {D^^)~^ (k) ~ (fe^)~°'^*^. 
The dominant term on the right-hand side is the ghost loop that originates from the action 
— TrlnM(S), whereas the subdominant terms — namely the tree term, the gluon loop and the 
two-loop term — all originate from the Yang- Mills action S'ym(S). 



D^uik) = (S^uk'^ -kt^kv) + (gluon loops) 

+ iV/(27r)"'^ / d'^p G{p + k) 




(8.2) 
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and tgh{x,y; B). The infrared asymptotic gluon and ghost propagators are designated 
{V)-^{B) = glldBy i^)xyiB) = fgh,a;y(5). The functional DS equation (6.10) for 
the ghost propagator is unchanged in form, as represented in Fig. 2, 

f y; B) ={-dH-' - gor'B^^d^) 5{x - y) 

, f , 6r2{u,y;B) (8-3) 

+ 9or'"'J dzdu Vil{x,z;B)d^g'^f{x,u;B)^^^^. 

In the infrared asymptotic hmit, only the ghost loop contributes to the functional DS 
equation for the gluon propagator (6.16), which reads 

iV-XU^,y;B) = -9or''{ J dzdu d^g'''ix,u;B) G-'^{x,z;B) (8.4) 

and is diagrammed in Fig. 4. 

An enormous simplification is apparent here, because the last equation allows an exact 
elimination of the asymptotic functional gluon propagator 'D'j^^{x, y; B). The one remaining 
unknown is the inverse ghost propagator rgh,xy(-B). 

When supplemented by the horizon condition (7.9), eqs. (8.3) and (8.4) are a finite 
system when expressed in terms of renormalized quantities. Indeed, with Zi = 1, renor- 
malization of the exact functional asymptotic equations is accomplished by writing 

B = Z^/^Br; t{B) = TRiBn); tg^{x, y- B) = Z-^Vg^^R{x, y- Br); 

go = {ZsZl^'^y^gR-, g{B) = Z^gR^BR); V^iu^x, y; B) = Z^^Vr^^^^x, y; Br). 

(8.5) 

Upon making these substitutions, the functional equation for the ghost propagator reads, 
r^lfl(^, y; Br) = {-dH-'z, - gRr'^'B^^^d^) 6{x - y) 

+ gnr"' J dzdu ^*;^,(a;, z; BR)d^g'^^{x, u; ^i^) ^^^^f^^^^^ > 

(8.6) 

where Z3 is given in (7.9), and the renormalized infrared asymptotic functional gluon 
propagator is given by 

/ /" - ^^i%h{u,z;BR)^tr 

{T^RXli^^y^BR) = -gRr'^ J dzdu d^g^^{x,u;BR) S%%x, z; BR)^g^^j^) . 

(8.7) 
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When (8.6) is expanded in a functional power series in B, appears only in the equation of 
order (5)° that determines the ghost propagator with source B = 0. This equation is finite 
as in the preceding section. All the higher order equations are independent of and finite. 
Equations (8.6) and (8.7) are a complete system of functional DS equations, diagrammed 
in Figs. 4 and 2, that are free of divergences, and that define the asymptotic infrared 
theory. The gluon propagator may be eliminated exactly from (8.7), and the asymptotic 
infrared theory is defined by the functional inverse ghost propagator Tgh,R{x, y; Bji). 

With Z^ given in (7.9), these equations are invariant under the finite renormalization- 
group transformations 

Br = zI'''b'r- Vr{Br) = V'r{B'r)- f gh,fl(a:, y; Br) = ~z^%^,r{x, y; Br)- 

QR = {z^zI'^^Qr- Qr{Br) = hQ'R{B'R)- V^^{x, y; Br) = Z3V'r^^^{x, y; B'r). 

(8.8) 

The quantity g\DR{k)G\{k) = gQD{k)G'^{k) is invariant under the renormalization (7.6) 
and (7.7). Consequently a scheme-independent running coupling constant, characteristic 
of the Landau gauge, may be defined [17] by Q;iand(^) = {4:'^)~^9oD{k)G'^{k){k'^)'^ . The 
asymptotic infrared theory is characterized, in addition to the infrared critical exponents 
ao and an, by aiand(O) ^ 8.915/A^, for color SU(N) [20]. 

The limit, in which the Yang- Mills action Sym{B) is systematically neglected, is a 
continuum analog of the lattice strong-coupling limit. Indeed if one rescales the gauge 
connection by the change of variable A' = QqA, the effective action, from which the DS 
equations were derived, reads 

E{A) = -TrlnM(A) + ^ym(A) 

= -TvlnM'iA') + {gl)-'Si,UA'), ^^'^^ 

where M'{A') and .SymI^') independent of go. Neglect of 5'ym(^) is the same as setting 
gQ^ = or, after renormalization, 51^^ = 0. The asymptotic infrared limit is described by 
the effective action 

E=-TrlnM(A). (8.10) 

A scheme- independent running coupling constant may be defined in the Coulomb gauge [69] 
by, q:coui(A;) = (47r)~^[12iV/(lliV - 2Nf)]k'^V(k), with Nf quarks in the fundamental representa- 
tion, where F(|fc|) = g'o li™fe4->oo -D44(fc, ^4). By contrast with aiand(fc) that is finite at fc = 0, it 
appears that acoui(fc) diverges like at small k, in a realization of infrared slavery that features 
a string tension, V{r) ~ cTcouif at large r [70] and [71]. 
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If one extends the non-perturbative formulation to a BRST-invariant theory, as out- 
hned in Appendix B, the BRST-invariant focal action reads 

S = Jd^x [s{df,cA^) + SYuiA)], (8.11) 

where the BRST operator acts according to 

sA^ = D^c; sc= sc=X; sA = 0. (8.12) 
The asymptotic infrared limit is described by the local BRST-invariant action 

S = Jd^x s{df,cA^) = j d^x (- d^c D^c + d^X A^), (8.13) 
and the infrared asymptotic correlators satisfy the Slavnov- Taylor identities. 

9. Mass gap 

The action S that describes the infrared asymptotic theory is not only BRST-invariant, 
it is BRST-exact, S = sX, and defines a topological quantum field theory. To see what its 
properties may be, recall that S describes the asymptotic infrared limit, in which external 
momenta k were small compared to Aqcdj so it is the limit Aqcd — > oo. If QCD is a 
theory with a mass gap of order Aqcd, then physical correlation lengths should vanish in 
the asymptotic theory, R ~ Aq^j-, — > 0. 

To show this, consider a gauge-invariant correlator, for example 

C{x) = (F^iA) Fi{A)) 

r (9 1) 

= N dAdcdcdX F^{A) F^{A) exp(-^), 
Jq. 

with a: 7^ 0, where Lorentz indices are suppressed F'^{x) F^y^{x)F'^j^{x), and the con- 
nected part is understood. Since the action is topological, we may make any transformation 
that commutes with s, without changing expectation values. As an example, consider the 
change of variable corresponding to a coordinate transformation — x'i^{x) oi A and c, 
leaving c and A unchanged. 



A'^{x') = ^^A,,{x); d{x') = c{x); c'{x) = c{x); )^{x) = X{x). (9.2) 
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(The result is the same if c and A are also transformed.) The infinitesimal form of this 
change of variable, with x'^ = x^ — ^^{x), is given by 

Af,{x) ^ A'^ix) = A^{x) + SA^ix) = A^{x) + ^^dxA^{x) + d^^^Ax{x) 
c{x) c{x) = c{x) + 5c{x) = c + c{x) + ^^dxc{x) (9.3) 
c{x) — > c'{x) — c{x); X{x) — > X'{x) — X{x). 

Upon making this change of variable in the functional integral, we obtain 

C{x)^ N [ dAdcdcdX F^{A') F^{A') exp[-S{A',c',c,X)], (9.4) 

where A' = A + 6A, and 

S{A', c', c. A) = J d^x (- a^c D^{A')c' + d^X A'^). (9.5) 

The integration in ^-space is cut-off at the Gribov horizon dQ' corresponding to M{A'). 
Integration over the ghost fields gives det M{A') which vanishes on the boundary dfl'. One 
may change the cut-off to the Gribov horizon dCt corresponding to M{A) because the error 
is only of order Moreover F^{A')Fq{A') is the coordinate transform of F^{A)Fq{A), 
which we write as 

f!{A') F^{A') = [1 + + m)]F!{A) F^{A), (9.6) 

where L{d^) is a numerical matrix that is linear in dx^^ and acts on the tensorial indices 
of F^{A) F§{A), and we have 

C{x) = [1 + ^^dx + L{d^)] N [ dAdcdcdX F^{A) F^{A) exp[-^(y4', c', c. A)]. (9.7) 

One may verify that s-operator commutes with the coordinate transformation, sA' = 
D{A')c', so 

S{A',c',c,X) = J d'^x s{d^cA'^) = S{A,c,c,X) + s6X, (9.8) 
where 6X = J d'^x d^c 5A^. Thus the variation of S is also s-exact, and we have 

C{x) = [1 + ^^dx + L{dC)] N I dAdcdcdX Fl{A) F^{A) (1 - s5X) exp[-^(A, c, c. A)] 

Jo. 

= [1 + eOx + L{dO\ { Fl{A) Fi{A) (1 - s5X) ). 

(9.9) 
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Gauge-invariant operators are s-invariant, 

( F^iA) Fi{A) s5X) = { s[F^{A) F^{A) 5X] ) = 0, (9.10) 

which vanishes because it is the expectation-value of an s-exact observable. This gives 
C{x) = [1 + i^d\ + L(9^)]C(x), so C{x) is invariant under arbitrary coordinate transfor- 
mation. Thus it is a number independent of x. It vanishes for x = oo, soG(a;) = for 
a; 0. The argument holds for a generic gauge-invariant correlator. 

We have shown that the correlation length R of gauge-invariant observables vanishes 
in the gauge- invariant, physical sector of asymptotic theory defined by S. In other words, 
the mass gap is infinite, M = 1/R = oo, in the physical sector of the asymptotic theory. It 
is tempting to conclude from this that there is a finite mass gap in the physical sector of the 
exact non-asymptotic theory, for otherwise we would have obtained non-zero correlators 
in the infrared limit. However local gauge-invariant observables like F'^{x) are composite 
operators, and so far we have discussed only the correlators of elementary fields. To 
establish that the mass gap in the non-asymptotic theory is finite, one should check that 
the correlators of local gauge-invariant operators in the limit of large separation are also 
given by the infrared asymptotic theory defined by S. 

10. Quarks 

So far we have neglected quarks, but they may be included in the time-independent 
Fokker-Planck equation [3]. The derivation of the non-perturbative Faddeev-Popov for- 
mula, including quarks, proceeds as in sees. 2 - 4, by changing quark variables according 
to '0 = g~^^ and ip = ^g. The result is that the quark action 5'qu = / d'^x '^{^^D^ +M)'0 
gets added to the gluon action E or S. According to the latest DS calculations that include 
Nf = 3 flavors of dynamical quarks, the quark-loop term in the DS-equation for gluons 
is subdominant in the infrared [72]. Provided that the effects of truncation are not too 
drastic, the quark contribution will also be subdominant in the infrared limit of the exact 
functional DS equation for the gluon propagator. In this case the inclusion of quarks does 
not disturb the simplicity of the gluon sector described by S. 

If the intrinsic mass of the quarks is finite, then the quark sector does not appear 
in the asymptotic infrared limit. If the instrinsic mass of the quarks is zero, the pion 
is a massless Goldstone boson associated with spontaneous breaking of chiral symmetry. 
However even in this case, in the (truncated) DS equation for the quark propagator given 
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in [72], the infrared limit of the quark propagator does not decouple from the degrees 
of freedom associated with finite momentum (in contrast to the gluon). This is to be 
expected because the parameters that characterize the dynamics of massless quarks, 
and /tt, are finite multiples of Aqcd, but the infrared asymptotic limit corresponds to 
Aqcd oo. Nevertheless one may ask if chiral symmetry is broken in the asymptotic 
infrared theory. The chiral-symmetry breaking parameter is given by {t/jt/j) = n{p{0, A)) , 
where p(A, A) is density of eigenvalues A, per unit volume, of the Dirac operator ■ D{A) 
in the configuration A. In the infrared asymptotic limit, the expectation-value (p(A, A)) is 
evaluated in the theory defined by the action S. One would expect that it gives ('0V') = OO; 
since this corresponds to Aqcd = oo. Thus in the theory defined by S, the average density 
of levels per unit volume (p(0,^)) of the Dirac operator • D{A) should be infinite at 
A = 0. 

The infrared asymptotic theory is far simpler than full QCD and provides a valu- 
able model in which the characteristic features of the confining phase, as described in 
the Landau gauge, are revealed. To understand confinement in the asymptotic theory, 
note that while the infrared components of A{x) are severely suppressed by the cut-off at 
the Gribov horizon, its short-wavc-lcngth components fluctuate wildly because the factor 
exp[— 5'ym(^)] is replaced by 1. Indeed, the infrared asymptotic gluon propagator D^^{k), 
eq. (8.1), is strongly enhanced in the ultraviolet. This suggests a picture of confinement 
in the infrared asymptotic theory in which the short-wave-length fluctuations of A°'{x) 
in color directions cause the decoherence of any field that carries a color charge. Indeed 
transport of a color vector q{T) along a path Zfj,{T), is described by Pexp(^o / AnZ^dt). 
In a highly random field ^^(a;), superposition of different paths is incoherent, so a field 
that bears a color charge does not propagate. In full QCD in Landau gauge, the dominant 
fiucuations of A{x) responsible for confinement should be on the length scale Aq^q. This 
picture of confinement is quite different from the scenario in Coulomb gauge, where confine- 
ment of color charge is attributed to a realization of infrared slavery by an instantaneous, 
long-range color-Coulomb potential [73], [70] and [71]. 

11. Conclusion 

We briefiy review the salient features of the non-perturbative continuum Euclidean 
formulation of QCD developed here. 
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(i) In Landau gauge one may integrate the Faddeev-Popov weight over the Gribov 
region Q instead of over the fundamental modular region A. 

(ii) The form of the Dyson-Schwinger equations is unchanged by the cut-off of the 
functional integral on the boundary dQ of the Gribov region, because the Faddeev-Popov 
determinent vanishes there. This simplicity makes the DS equations the method of choice 
for non-perturbative calculations in QCD. 

(iii) The restriction to the Gribov region provides supplementary conditions that gov- 
ern the choice of solution of the DS equations. Two conditions are the positivity of the gluon 
and ghost propagators. Another is the horizon condition which is the statement that the 
ghost propagator G{k) is more singular than in the infrared, limk->o[k'^G{k)]~^ = 0. 
This fixes the ghost-propagator renormalization constant Z3 to the value (7.9). Al- 
though (7.9) is in flagrant disagreement with the perturbative expression for Z3, nev- 
ertheless it is consistent with the perturbative renormalization group. 

(iv) Implementation of the horizon condition in the DS equations puts QCD into a 
non-perturbative phase. 

(v) Recent solutions of the truncated DS equations possess an asymptotic infrared 
limit that is obtained by systematically neglecting the terms in the DS equations that 
come from the Yang-Mills action S'ym(^), but keeping the Faddeev-Popov determinant 
and the cut-off at the Gribov horizon. If the effects of truncation are not too drastic, this 
also gives an exact asymptotic infrared limit of QCD that is a continuum analog of the 
strong-coupling limit in lattice gauge theory. This is possible because convergence of the 
A-integration without the Yang-Mills factor exp[— 5'ym(^)] may be assured by the cut-off 
at the Gribov horizon. 

(vi) The asymptotic infrared limit of QCD is defined by the functional DS equations 
(8.6) and (8.7). The gluon propagator may be eliminated exactly from (8.7), and the 
asymptotic infrared theory is completely characterized by the functional inverse ghost 
propagator Tgh{x, y',B). 

(vii) There exists a local BRST-invariant extension of the present non-perturbative for- 
mulation, sketched out in Appendix B. This ensures that the Slavnov- Taylor identities hold 
in the non-perturbative theory. The asymptotic infrared limit of QCD, valid at distances 
large compared to 1/Aqcd, is described by the BRST-exact action, S = J d^x s {d^cA,j)^ 
that defines a topological quantum field theory with an infinite mass gap. 

36 



(viii) The extension of the non-perturbative formulation to include the quark action 
J d^x 'ip{'yij,D^+m)'ip is immediate. The presence of quarks does not disturb the asymptotic 
infrared limit of the gluon sector. 

(ix) The asymptotic infrared theory provides a simple model in the Landau gauge 
in which the characteristic features of confinement may be understood. A picture of 
confinement of color charge emerges, in which the highly random fluctuations of the gluon 
field A cause the superposition from the transport of color charge along different paths to 
interfere incoherently, so the fields that bear a color charge do not propagate. 

Acknowledgment s 

It is a pleasure to thank Reinhard Alkofer, Laurent Baulieu, Alexander Rutenburg, 
S. R. S. Varadhan, and Lorenz von Smekal for valuable discussions and correspondence. 
This research was partially supported by the National Science Foundation under grant 
PHY-0099393. 



Appendix A. Resolution of paradox 

At first sight it is surprising that expectation-values taken over the fundamental mod- 
ular region A and the Gribov region Q are equal. In this Appendix we show how this 
paradox is resolved. 

A.l. Argument of Semenov-Tyan-Shanskii and Pranke 

The proof by Semenov-Tyan-Shanskii and Franke [4] that the Gribov region Q and 
the fundamental modular region A are different, substantiated by instances are given 
in [6], was long considered to disprove (1.3). We review the argument of [4]. Let 
g{t) = exp{tu}) be a one-parameter subgroup of the local gauge group with generator 
Lu = uj{x). To be definite, we normalize uj to (u^lo) = V, where V is the Euclidean 
volume. Let Af^{t,uj, B) = g{t)~^B^g{t) -\- g{t)~^d^g{t), be the gauge-transform of 
under g{t) = exp(tc<;), so A{0,uj,B) = B, and let FB{t,Lv) be the Hilbert square norm of 
A{t, u>, B), regarded as a function of t and u> for fixed B, 

FB{t,uj) = \\A{t,uj,B)\\'' = J d'^x \A^{t,u;,B)\\ (A.l) 
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The fundamental modular region A is the set of B such that ^^(0,0;) is an absolute 
minimum, Fb{0,(jj) < FB{t,uj) for all oj and t. The Gribov region Cl is the set of B for 
which Fb(0, a;) is a relative minimum Fb(0, a;) < Fsit, cj) for all uj and sufficiently small t. 
We differentiate FB{t,oj) with respect to t, and use A'^ = Dfj,{A)uj = -D^a;, 

F's{t,u;) = 2 {D^uj,A^) = 2{d^u;,A^) = -2{u;,d^A^) 
F'^{t,uj) = 2 (a^a;,D^a;) = -2 {u.d^D^u) 

(A.2) 

F^'(t,a;) = 2 (^^o;, D^a; x a;) 

cc;) = 2 (a^c.;, (I^^cc; x c^) x u), 

where X = [X, Y] is the commutator in the Lie algebra. These formulas show that the 
interior of VL consists of all transverse configurations B, d ■ B = 0, such that all non-trivial 
eigenvalues of M{B) = —d^D^{B) are strictly positive, Xn{B) > 0. Moreover for B on 
the boundary dO,^ M{B) has at least one non-trivial eigenvalue that vanishes, Xi{B) = 0. 

We specialize to the SU(2) group, so the commutator X x y is the ordinary 3-vector 
cross product. The vector triple product gives 

F'^'it, uj) =2 (df^u, {uj ■ D^iv ix) -ij? Df^uj) ) 

= 2 {d^io, CO ■ d^u u)) +2 (w, d^{u'^D^u)) 

= 2{u ■ d^cu, UJ ■ d^cj) + 2{uj, d^{uj^)D^uj) + 2{uj, oj^d^D^u) (A.3) 
= (l/2)(a^(a;2), d^iuj^)) + 2(a; • d^uj, d^iu^)) + 2{uj\ uj ■ d^D^uj) 

= {3/2){d^{uj^),d^{u^)) + 2{u\ UJ ■ d^D^uj), 

where the dot is contraction on color indices. 

Let i? be a point on the Gribov horizon dfl, so B is transverse df^Bjj^ = 0, and the 
Faddeev-Popov operator —dy^Dn^B) is non- negative, but with at least one non-trivial null 
eigenvalue, d^D^{B)uo = 0, for some ujq. By (A.2), we have Fg(0,a;o) = F'^{0,ujo) = 0. 
For B on dQ, it follows that in general Fb{0, uj) is not a local minimum on the gauge orbit 
through B because, in general, F'^'{0,ujo) ^ 0, so FB{t,uJo) — Fb{0,ujo) changes sign at 
t = 0. By continuity this implies that nearby points inside the Gribov region Q cannot be 
absolute minima, even though they are relative minima. They are Gribov copies inside fl. 
This is the argument of [4], and examples for which F^'(0, cuq) ^ 0, are given in [6]. 

But let's evaluate the 4th derivative at t = 0, in the direction ujq. With d^D,j_{B)ujQ = 
0, we have from (A.3), 

Ff'(0,a;o) = (3/2) J d^x [d^{uj'o)]'. (A.4) 
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This is the integral of a positive density, and we expect that F^"(0, a;o) is large and positive. 

The relevant question for comparing the expectation values over Q and over A is not 
whether these regions coincide — they do not — but whether the normalized averages 
over these sets are equal in the thermodynamic limit. Here we implicitly suppose a lattice 
discretization, and configurations that are sampled from the Wilson ensemble. In the 
thermodynamic limit, the probability may get concentrated on a subset that consists of 
a boundary or part of a boundary. The boundaries of A and O may approach each other 
in the thermodynamic limit for typical configurations on the boundary. If F'^"{0,ujo) is 
large, and F'^'{0,u}o) is small, then there is a local minimum near B, which could be the 
absolute minimum on the gauge orbit. If the distance to the absolute minimum vanishes 
in the thermodynamic limit for a typical configuration, then the argument of [4] does not 
disprove (1.3). 

We normalize cuq to (a;o,a;o) = V, where V is the volume of Euclidean space. We 
estimate quantities using this normalization, and we shall verify that the conclusions do 
not depend on the normalization of coq. With this normalization, we estimate that u>o{x) = 
0(1). Since F^"(0,c<;o) is the integral of a positive local density over a volume V, we 
estimate that F'^"{0, ujq) = 0{V), for a typical configuration B on the Gribov horizon. On 
the other hand, the density that appears in F^'(0, cjq) has no definite sign. For a typical 
configuration, sampled from the Wilson ensemble, we make the crudest statistical estimate 
namely random density, so F'^'(0,ujq) = 0(F^/^). This is small compared to F^"(0,a;o). 
We seek a nearby minimum on the gauge orbit through B. For simplicity we assume that 
all non-trivial eigenvalues of M{B) are strictly positive, apart from the zero eigenvalue 
belonging to uq, which is the only dangerous direction. We write F{t) = FB{t,oJo), and 
we have 

F{t) = F(0) + (l/3!)F'"(0)t3 + (l/4!)F""(0)^^ (A.5) 

with neglect of higher order terms. The minimum is found at F'{tci) = 0, which gives 
t„ = -3F'"(0)/F""(0), and one has, 

F(U = F(0)-(9/8)i|;^. (A.6) 

This is lower than F{0), in agreement with the argument of [4]. This expression is indepen- 
dent of the normalization of uo, as one sees from (A. 2), so our estimate for this quantity 
is independent of the normalization of uq. By the above estimates, the second term is of 
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order = V ^. It is small compared to the first term, F(0) = \ \B\\^, which is 

of order V. The configuration at the nearby minimum is 

3F"'(0) (A.7) 
= ^^(^)- [D,iB)uo]{x), 

which is again independent of the normalization of uq. According to the above estimates, 
the second term is of order V~^/'^. Thus in the thermodynamic limit of lattice gauge 
theory, V ^ oo, the nearby minimum approaches the point B on the Gribov horizon. 
In actuality, the problem of minimizing the functional F^ig) = on the lattice is 

a problem of spin-glass type, so one expects many, nearly degenerate, relative minima, 
and the one found here is not necessarily the absolute minimum. Nevertheless the point 
remains that [4] does not disprove the equality of expectation-values on A and Q in the 
thermodynamic limit. 

A. 2. Many Gribov copies inside the Gribov region from numerical simulations 

We now consider the fact that in numerical gauge-fixing to Landau gauge in lattice 
gauge theory, there are many local minima (i.e. Gribov copies inside the Gribov region, Q), 
on a typical gauge orbit, [7], [8], [9], [10], [11]. Their number grows with the lattice size as 
is characteristic of a spin-glass. In this sense is very large compared to A. However the 
number of dimensions of configuration space is high, and our geometrical intuition from 
3-space may be misleading. Indeed, on a lattice of Euclidean volume V, the dimension D 
of configuration space is D = fV, where / is the number of degrees of freedom per lattice 
site, and the dimension D of configuration space diverges with the Euclidean volume V. 

In continuum gauge theory A and O are both convex and bounded in every direc- 
tion [4]. By simple entropy considerations, the population in a bounded region of a high- 
dimensional space gets concentrated on the boundary. For example inside a sphere of 
radius i? in a D-dimensional space, the radial density is given by r^~^dr, and for r < i? is 
highly concentrated near the boundary r = R. To take the simplest example, consider two 
spheres (in configuration space), the first of radius -R, and the second of radius R + cV~^/'^. 
In the spirit of the previous estimates, these would be the radii of A and Q. The ratio of 
the radii {R + cV~^/'^)/R approaches unity, in the limit F — > oo, so all n-th moments, (r"') 
for finite n, of the two spheres become equal. On the other hand the ratio of their volumes 
is given by [{R + cV-'^/^)/R]^ = [{R+ cV-^/^)/R]f^ , where D = fV is the dimension 
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of configuration space. For large V the ratio of the volumes of the two spheres is thus 
exp(a/F^/^/i?), which diverges exponentially like V^^'^. In this example the ratio of the 
volumes of the two spheres diverges with V, but all finite moments of the two spheres be- 
come equal! In field theory the n-th moments of the distribution are the n-point functions 
{A{xi)...A{xn)- So again, the fact that there are many Gribov copies inside il, does not 
disprove that averages calculated over A or f2 are equal. 

A. 3. Gauge theory on a finite lattice 

For a finite lattice the paradox becomes acute. Stochastic quantization may also be 
defined in lattice gauge theory [10]. As in the continuum theory, a drift force a~^Kgi^ 
tangent to the gauge orbit may be chosen in the direction of steepest descent of a suitable 
minimizing function, and is globally resoring. It appears that one may solve the lattice 
Fokker-Planck equation in the limit a — > on a finite lattice, by the method used in sees. 2 
to 4, for it depends only on general geometrical properties that are common to lattice and 
continuum gauge theories. If so, one would again be led to the conclusion that the weight 
inside the Gribov region is given by the lattice analog of (4.17) namely ATexpf— S'w(C/)], 
where Sw{U) is the Wilson action, and U is a configuration in the lattice Gribov region Q. 
However on a finite lattice the distinction between the fundamental modular region A and 
the Gribov region Q can surely not be ignored. The resolution of this paradox would 
appear to be that in lattice gauge theory the Gribov region ft is made of disconnected 
pieces ilj. In each piece, the solution is indeed given by Qi{U) = Nie-xp[—Sw{U)], for 
U E fli, where the normalizations Ni are left indeterminate by the method of sees. 2 to 4. 
Presumably, the average with the lattice Faddeev-Popov weight over all the disconnected 
pieces Qi of the Gribov region, with the correct normalization Ni in each piece, will agree 
with with same integral over the fundamental modular region A. 

Appendix B. BRST-invariant formulation 

New issues arise when the non-perturbative approach is extended to a theory with a 
local BRST-invariant action. 
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B.l. Off-shell transversality condition 

To obtain a local action, on must take the transversality condition "off-shell". The 
off-shell partition function is given by 

Z{J,L)= I DADX detM{A) exp[-5YM(^) + i(A, 9 • A) + i( J, A) + i(L, A)], (B.l) 
Jo. 

where A is the Nakanishi-Lautrup Lagrange multiplier field that enforces the gauge con- 
dition 5-^ = 0, and L is its source. This reduces to (5.1) for L = 0. It is not immedi- 
ately obvious what region Q, to integrate over because A is not transverse for L 7^ 0, so 
M{A) = —d ■ D{A) is not a symmetric operator. One must also take the Gribov hori- 
zon dO, off shell when the gauge condition is off-shell. If we effect the A integration, the 
last integral becomes, 

Z(J,L)= / DA detM{A) 6{d-A + L) exp[-SYM{A) + i{J,A)]. (B.2) 

Only configurations A of the form A = B — d{d'^)~^L are relevant, where B is trans- 
verse. We regard the partition function Z( J, L) as a formal power series in the source L. 
Both the lowest non-trivial eigenvalue \\[B — d{d'^)~^L] of the Faddeev-Popov operator, 
M[B — 5(5^)~^L], and the points Bq{L) where it vanishes, may be calculated by formal 
perturbation theory as a power series in L. Here -Bo(O) is a point on the on-shell horizon. 
In this way we may take the Gribov horizon dO. off-shell. 

B.2. Faddeev-Popov ghosts 

One may make the action local by writing 

detM(S) = J DcDc e}cp{c,M{B)c), (B.3) 

where c and c are anti-commuting ghost and anti-ghost fields. Grassmannian sources, 77 
and f], are then introduced, so this gets replaced by 

J DcDc exp[(c,M(S)c) + (f7,c) + {c,r])] = det M{B)exp{fj,M-^{B)r]). (B.4) 

This expression does not vanish on the boundary dfl. For by an eigenfunction expansion 
of M~^{B), we obtain for the last expression 

n exp [ ^?7n?7n] = H 11 ("^ + ^^n?7n) = JI ('^^ + ^"^'^) " (B-5) 
n ^ n n ^ 
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It does not contain as a factor Xi{B) that vanishes on dfl. For this reason, we did not use 
Faddeev-Popov ghost fields and their sources in the derivation of the DS equations in sees. 5 
and 6. Nevertheless we obtained the same DS equations, including the ghost propagators, 
that we would have obtained if we had introduced the ghost fields and their sources. For 
this reason, and by use of the off-shell Gribov horizon, it should be possible to extend 
the non-perturbative approach to the theory defined by the familiar BRST-invariant local 
action (8.11), integrated over the off-shell Gribov region. 

Appendix C. Properties of the Gribov region 

We note three properties of the Gribov region Q defined in (1.2). (i) O contains the 
origin A = 0. (ii) It is bounded in every direction, (iii) It is convex. We give the one-line 
proofs of these properties [74]. They follow from the expression M[A) = Mq + Mi{A), 
where Mq'^{A) = — 9^5""^, and Mf '^(A) = — gfo/^^^A^c?^, where A is transverse. Property (i) 
is obvious since Mq = is strictly positive. To establish (ii), note that Mi{A) has 

zero trace, since it is traceless on color indices = 0. Thus, for any given ^, there exists 
a state uj for which the expectation value of Mi{A) is is negative, E = {uj, Mi{A)(jj) < 0. 
Moreover Mi{A) is linear in A, Mi{XA) = XMi{A), so upon replacing A by A^, where A is 
a positive number, we have (a;, M{XA)(jj) = (a;, Moa;)-|-A(a;, Mi{A)(jj) = (a;, Moa;)-|-A£'. By 
taking A sufficiently large and positive, the expectation value is negative (a;, M{XA)(jj) < 0. 
This establishes (ii). To establish convexity, we must show that M{aAi+PA2) is a strictly 
positive operator when M{Ai) and M{A2) are both strictly positive operators, for all 
positive a and /3, with a + (3 = 1. This is immediate because Mi{A) depends linearly on 
A, and we have M{aAi + (3A2) = aM{Ai) + l3M{A2). QED 
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Figure Captions 



Fig. 1. The coordinate patch U in ^-space is the clam-shaped region viewed edge on. 
The Gribov region Q is represented by the thick horizontal line. 

Fig. 2. The functional DS equation (6.10) for the complete ghost propagator Q{x,y;B) 
in the presence of the source B. The thin line is the tree-level term. The heavy 
line with (without) the arrow is the complete ghost (gluon) propagator Q{x, y; B) 
{V{x,y;B)) in the presence of the source B. The circle is the complete ghost- 
ghost-gluon vertex in the presence of the source B. 

Fig. 3. The functional DS equation (6.16) for the complete gluon propagator T>{x,y;B) 
in the presence of the source B. The thin line is the tree-level term. The heavy 
line with (without) the arrow is the complete ghost (gluon) propagator Q{x,y; B) 
{'D{x,y;B)) in the presence of the source B. The circles are complete 3- and 
4- vertices in the presence of the source B. 

Fig. 4. The functional DS equation (8.4) for the complete infrared asymptotic gluon 
propagator V{x^ y; B) in the presence of the source B. There is no tree term nor 
any gluon loop, but only the ghost loop. The heavy line with the arrow is the 
complete infrared asymptotic ghost propagator Q{x, y; B) in the presence of the 
source B. The functional DS equation (8.6) for the complete infrared asymptotic 
ghost propagator Q{x, y; B) in the presence of the source B is as in Fig. 2. 
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D(-1)(x,y;B)= x y 




